Binding and Activation of LRP1-Dependent Cell Signaling in Schwann Cells Using a Peptide Derived from the Hemopexin Domain of MMP-9.

Schwann cells (SCs) undergo phenotypic transformation and then orchestrate nerve repair following a peripheral nervous system injury. The low-density lipoprotein receptor-related protein-1 (LRP1) is significantly upregulated in SCs in response to acute injury, activating cJun and promoting SC survival. Matrix-metalloproteinase-9 (MMP-9) is an LRP1 ligand that binds LRP1 through its hemopexin domain (PEX) and activates SC survival signaling and migration. To identify novel peptide mimetics within the hemopexin domain of MMP-9, we examined the crystal structure of PEX, synthesized four peptides, and examined their potential to bind and activate LRP1. We demonstrate that a 22 amino acid peptide, peptide 2, was the only peptide that activated Akt and ERK1/2 signaling in SCs, similar to a glutathione s-transferase (GST)-fused holoprotein, GST-PEX. Intraneural injection of peptide 2, but not vehicle, into crush-injured sciatic nerves activated cJun greater than 2.5-fold in wild-type mice, supporting that peptide 2 can activate the SC repair signaling in vivo. Peptide 2 also bound to Fc-fusion proteins containing the ligand-binding motifs of LRP1, clusters of complement-like repeats (CCRII and CCRIV). Pulldown and computational studies of alanine mutants of peptide 2 showed that positively charged lysine and arginine amino acids within the peptide are critical for stability and binding to CCRII. Collectively, these studies demonstrate that a novel peptide derived from PEX can serve as an LRP1 agonist and possesses qualities previously associated with LRP1 binding and SC signaling in vitro and in vivo.

Integrating cellular electron microscopy with multimodal data to explore biology across space and time.

Biology spans a continuum of length and time scales. Individual experimental methods only glimpse discrete pieces of this spectrum but can be combined to construct a more holistic view. In this Review, we detail the latest advancements in volume electron microscopy (vEM) and cryo-electron tomography (cryo-ET), which together can visualize biological complexity across scales from the organization of cells in large tissues to the molecular details inside native cellular environments. In addition, we discuss emerging methodologies for integrating three-dimensional electron microscopy (3DEM) imaging with multimodal data, including fluorescence microscopy, mass spectrometry, single-particle analysis, and AI-based structure prediction. This multifaceted approach fills gaps in the biological continuum, providing functional context, spatial organization, molecular identity, and native interactions. We conclude with a perspective on incorporating diverse data into computational simulations that further bridge and extend length scales while integrating the dimension of time.

Glycoproteomic landscape and structural dynamics of TIM family immune checkpoints enabled by mucinase SmE.

Mucin-domain glycoproteins are densely O-glycosylated and play critical roles in a host of biological functions. In particular, the T cell immunoglobulin and mucin-domain containing family of proteins (TIM-1, -3, -4) decorate immune cells and act as key regulators in cellular immunity. However, their dense O-glycosylation remains enigmatic, primarily due to the challenges associated with studying mucin domains. Here, we demonstrate that the mucinase SmE has a unique ability to cleave at residues bearing very complex glycans. SmE enables improved mass spectrometric analysis of several mucins, including the entire TIM family. With this information in-hand, we perform molecular dynamics (MD) simulations of TIM-3 and -4 to understand how glycosylation affects structural features of these proteins. Finally, we use these models to investigate the functional relevance of glycosylation for TIM-3 function and ligand binding. Overall, we present a powerful workflow to better understand the detailed molecular structures and functions of the mucinome.

Multiscale computational modeling of the effects of 2'-deoxy-ATP on cardiac muscle calcium handling.

2'-Deoxy-ATP (dATP), a naturally occurring near analog of ATP, is a well-documented myosin activator that has been shown to increase contractile force, improve pump function, and enhance lusitropy in the heart. Calcium transients in cardiomyocytes with elevated levels of dATP show faster calcium decay compared with cardiomyocytes with basal levels of dATP, but the mechanisms behind this are unknown. Here, we design and utilize a multiscale computational modeling framework to test the hypothesis that dATP acts on the sarcoendoplasmic reticulum calcium-ATPase (SERCA) pump to accelerate calcium re-uptake into the sarcoplasmic reticulum during cardiac relaxation. Gaussian accelerated molecular dynamics simulations of human cardiac SERCA2A in the E1 apo, ATP-bound and dATP-bound states showed that dATP forms more stable contacts in the nucleotide binding pocket of SERCA and leads to increased closure of cytosolic domains. These structural changes ultimately lead to changes in calcium binding, which we assessed using Brownian dynamics simulations. We found that dATP increases calcium association rate constants to SERCA and that dATP binds to apo SERCA more rapidly than ATP. Using a compartmental ordinary differential equation model of human cardiomyocyte excitation-contraction coupling, we found that these increased association rate constants contributed to the accelerated rates of calcium transient decay observed experimentally. This study provides clear mechanistic evidence of enhancements in cardiac SERCA2A pump function due to interactions with dATP.

Revealing the Impacts of Chemical Complexity on Submicrometer Sea Spray Aerosol Morphology.

Sea spray aerosol (SSA) ejected through bursting bubbles at the ocean surface is a complex mixture of salts and organic species. Submicrometer SSA particles have long atmospheric lifetimes and play a critical role in the climate system. Composition impacts their ability to form marine clouds, yet their cloud-forming potential is difficult to study due to their small size. Here, we use large-scale molecular dynamics (MD) simulations as a "computational microscope" to provide never-before-seen views of 40 nm model aerosol particles and their molecular morphologies. We investigate how increasing chemical complexity impacts the distribution of organic material throughout individual particles for a range of organic constituents with varying chemical properties. Our simulations show that common organic marine surfactants readily partition between both the surface and interior of the aerosol, indicating that nascent SSA may be more heterogeneous than traditional morphological models suggest. We support our computational observations of SSA surface heterogeneity with Brewster angle microscopy on model interfaces. These observations indicate that increased chemical complexity in submicrometer SSA leads to a reduced surface coverage by marine organics, which may facilitate water uptake in the atmosphere. Our work thus establishes large-scale MD simulations as a novel technique for interrogating aerosols at the single-particle level.

Selectivity and Ranking of Tight-Binding JAK-STAT Inhibitors Using Markovian Milestoning with Voronoi Tessellations.

Janus kinases (JAK), a group of proteins in the nonreceptor tyrosine kinase (NRTKs) family, play a crucial role in growth, survival, and angiogenesis. They are activated by cytokines through the Janus kinase-signal transducer and activator of a transcription (JAK-STAT) signaling pathway. JAK-STAT signaling pathways have significant roles in the regulation of cell division, apoptosis, and immunity. Identification of the V617F mutation in the Janus homology 2 (JH2) domain of JAK2 leading to myeloproliferative disorders has stimulated great interest in the drug discovery community to develop JAK2-specific inhibitors. However, such inhibitors should be selective toward JAK2 over other JAKs and display an extended residence time. Recently, novel JAK2/STAT5 axis inhibitors (N-(1H-pyrazol-3-yl)pyrimidin-2-amino derivatives) have displayed extended residence times (hours or longer) on target and adequate selectivity excluding JAK3. To facilitate a deeper understanding of the kinase-inhibitor interactions and advance the development of such inhibitors, we utilize a multiscale Markovian milestoning with Voronoi tessellations (MMVT) approach within the Simulation-Enabled Estimation of Kinetic Rates v.2 (SEEKR2) program to rank order these inhibitors based on their kinetic properties and further explain the selectivity of JAK2 inhibitors over JAK3. Our approach investigates the kinetic and thermodynamic properties of JAK-inhibitor complexes in a user-friendly, fast, efficient, and accurate manner compared to other brute force and hybrid-enhanced sampling approaches.

Inhibition of PRMT5 by market drugs as a novel cancer therapeutic avenue.

Market drugs, such as Food and Drug Administration (FDA) or European Medicines Agency (EMA)-approved drugs for specific indications provide opportunities for repurposing for newer therapeutics. This potentially saves resources invested in clinical trials that verify drug safety and tolerance in humans prior to alternative indication approval. Protein arginine methyltransferase 5 (PRMT5) overexpression has been linked to promoting the tumor phenotype in several cancers, including pancreatic ductal adenocarcinoma (PDAC), colorectal cancer (CRC), and breast cancer (BC), making PRMT5 an important target for cancer therapy. Previously, we showed that PRMT5-mediated methylation of the nuclear factor (NF)-κB, partially contributes to its constitutive activation observed in cancers. In this study, we utilized an AlphaLISA-based high-throughput screening method adapted in our lab, and identified one FDA-approved drug, Candesartan cilexetil (Can, used in hypertension treatment) and one EMA-approved drug, Cloperastine hydrochloride (Clo, used in cough treatment) that had significant PRMT5-inhibitory activity, and their anti-tumor properties were validated using cancer phenotypic assays in vitro. Furthermore, PRMT5 selective inhibition of methyltransferase activity was confirmed by reduction of both NF-κB methylation and its subsequent activation upon drug treatment. Using in silico prediction, we identified critical residues on PRMT5 targeted by these drugs that may interfere with its enzymatic activity. Finally, Clo and Can treatment have exhibited marked reduction in tumor growth in vivo. Overall, we provide basis for pursuing repurposing Clo and Can as anti-PRMT5 cancer therapies. Our study offers potential safe and fast repurposing of previously unknown PRMT5 inhibitors into clinical practice.

SARS-CoV-2 evolved variants optimize binding to cellular glycocalyx.

Viral variants of concern continue to arise for SARS-CoV-2, potentially impacting both methods for detection and mechanisms of action. Here, we investigate the effect of an evolving spike positive charge in SARS-CoV-2 variants and subsequent interactions with heparan sulfate and the angiotensin converting enzyme 2 (ACE2) in the glycocalyx. We show that the positively charged Omicron variant evolved enhanced binding rates to the negatively charged glycocalyx. Moreover, we discover that while the Omicron spike-ACE2 affinity is comparable to that of the Delta variant, the Omicron spike interactions with heparan sulfate are significantly enhanced, giving rise to a ternary complex of spike-heparan sulfate-ACE2 with a large proportion of double-bound and triple-bound ACE2. Our findings suggest that SARS-CoV-2 variants evolve to be more dependent on heparan sulfate in viral attachment and infection. This discovery enables us to engineer a second-generation lateral-flow test strip that harnesses both heparin and ACE2 to reliably detect all variants of concern, including Omicron.

APOBEC3A Catalytic Inactivity Mutation Induces Tertiary Structure Destabilization.

APOBEC3A (A3A)-catalyzed DNA cytosine deamination is implicated in virus and cancer mutagenesis, and A3A is a target for small molecule drug discovery. The catalytic glutamic acid (E72) is frequently mutated in biochemical studies to characterize deamination-dependent versus deamination-independent A3A functions. Additionally, catalytically active A3A is toxic in bacterial expression systems, which adversely affects yield during recombinant A3A expression. Here, we demonstrate that mutating the catalytic glutamic acid to an isosteric glutamine (E72Q) significantly decreases the thermal stability of the protein, compared to the alanine-inactivating mutation (E72A). Differential scanning fluorimetry and mass spectrometry reveal that A3A E72Q is less thermally stable than A3A E72A or wild-type A3A. Strikingly, A3A E72Q is partially denatured at 37 °C and binds single-stranded DNA with significantly poorer affinity compared to A3A E72A. This study constitutes an important cautionary note on A3A protein design and informs that A3A E72A is the preferred catalytic inactivation mutation for most applications.

Transmissible SARS-CoV-2 variants with resistance to clinical protease inhibitors.

Vaccines and drugs have helped reduce disease severity and blunt the spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). However, ongoing virus transmission, continuous evolution, and increasing selective pressures have the potential to yield viral variants capable of resisting these interventions. Here, we investigate the susceptibility of natural variants of the main protease [Mpro; 3C-like protease (3CLpro)] of SARS-CoV-2 to protease inhibitors. Multiple single amino acid changes in Mpro confer resistance to nirmatrelvir (the active component of Paxlovid). An additional clinical-stage inhibitor, ensitrelvir (Xocova), shows a different resistance mutation profile. Importantly, phylogenetic analyses indicate that several of these resistant variants have pre-existed the introduction of these drugs into the human population and are capable of spreading. These results encourage the monitoring of resistance variants and the development of additional protease inhibitors and other antiviral drugs with different mechanisms of action and resistance profiles for combinatorial therapy.

Glycoproteomic landscape and structural dynamics of TIM family immune checkpoints enabled by mucinase SmE.

Mucin-domain glycoproteins are densely O-glycosylated and play critical roles in a host of biological functions. In particular, the T cell immunoglobulin and mucin-domain containing family of proteins (TIM-1, -3, -4) decorate immune cells and act as key checkpoint inhibitors in cancer. However, their dense O-glycosylation remains enigmatic both in terms of glycoproteomic landscape and structural dynamics, primarily due to the challenges associated with studying mucin domains. Here, we present a mucinase (SmE) and demonstrate its ability to selectively cleave along the mucin glycoprotein backbone, similar to others of its kind. Unlike other mucinases, though, SmE harbors the unique ability to cleave at residues bearing extremely complex glycans which enabled improved mass spectrometric analysis of several mucins, including the entire TIM family. With this information in-hand, we performed molecular dynamics (MD) simulations of TIM-3 and -4 to demonstrate how glycosylation affects structural features of these proteins. Overall, we present a powerful workflow to better understand the detailed molecular structures of the mucinome.

DeepWEST: Deep Learning of Kinetic Models with the Weighted Ensemble Simulation Toolkit for Enhanced Sampling.

Recent advances in computational power and algorithms have enabled molecular dynamics (MD) simulations to reach greater time scales. However, for observing conformational transitions associated with biomolecular processes, MD simulations still have limitations. Several enhanced sampling techniques seek to address this challenge, including the weighted ensemble (WE) method, which samples transitions between metastable states using many weighted trajectories to estimate kinetic rate constants. However, initial sampling of the potential energy surface has a significant impact on the performance of WE, i.e., convergence and efficiency. We therefore introduce deep-learned kinetic modeling approaches that extract statistically relevant information from short MD trajectories to provide a well-sampled initial state distribution for WE simulations. This hybrid approach overcomes any statistical bias to the system, as it runs short unbiased MD trajectories and identifies meaningful metastable states of the system. It is shown to provide a more refined free energy landscape closer to the steady state that could efficiently sample kinetic properties such as rate constants.

Targeted protein S-nitrosylation of ACE2 inhibits SARS-CoV-2 infection.

Prevention of infection and propagation of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a high priority in the Coronavirus Disease 2019 (COVID-19) pandemic. Here we describe S-nitrosylation of multiple proteins involved in SARS-CoV-2 infection, including angiotensin-converting enzyme 2 (ACE2), the receptor for viral entry. This reaction prevents binding of ACE2 to the SARS-CoV-2 spike protein, thereby inhibiting viral entry, infectivity and cytotoxicity. Aminoadamantane compounds also inhibit coronavirus ion channels formed by envelope (E) protein. Accordingly, we developed dual-mechanism aminoadamantane nitrate compounds that inhibit viral entry and, thus, the spread of infection by S-nitrosylating ACE2 via targeted delivery of the drug after E protein channel blockade. These non-toxic compounds are active in vitro and in vivo in the Syrian hamster COVID-19 model and, thus, provide a novel avenue to pursue therapy.

Spike-heparan sulfate interactions in SARS-CoV-2 infection.

Recent biochemical, biophysical, and genetic studies have shown that heparan sulfate, a major component of the cellular glycocalyx, participates in infection of SARS-CoV-2 by facilitating the so-called open conformation of the spike protein, which is required for binding to ACE2. This review highlights the involvement of heparan sulfate in the SARS-CoV-2 infection cycle and argues that there is a high degree of coordination between host cell heparan sulfate and asparagine-linked glycans on the spike in enabling ACE2 binding and subsequent infection. The discovery that spike protein binding and infection depends on both viral and host glycans provides insights into the evolution, spread and potential therapies for SARS-CoV-2 and its variants.

Discovery of compounds that reactivate p53 mutants in vitro and in vivo.

The tumor suppressor p53 is the most frequently mutated protein in human cancer. The majority of these mutations are missense mutations in the DNA binding domain of p53. Restoring p53 tumor suppressor function could have a major impact on the therapy for a wide range of cancers. Here we report a virtual screening approach that identified several small molecules with p53 reactivation activities. The UCI-LC0023 compound series was studied in detail and was shown to bind p53, induce a conformational change in mutant p53, restore the ability of p53 hotspot mutants to associate with chromatin, reestablish sequence-specific DNA binding of a p53 mutant in a reconstituted in vitro system, induce p53-dependent transcription programs, and prevent progression of tumors carrying mutant p53, but not p53null or p53WT alleles. Our study demonstrates feasibility of a computation-guided approach to identify small molecule corrector drugs for p53 hotspot mutations.

Architecture and self-assembly of the jumbo bacteriophage nuclear shell.

Bacteria encode myriad defences that target the genomes of infecting bacteriophage, including restriction-modification and CRISPR-Cas systems1. In response, one family of large bacteriophages uses a nucleus-like compartment to protect its replicating genomes by excluding host defence factors2-4. However, the principal composition and structure of this compartment remain unknown. Here we find that the bacteriophage nuclear shell assembles primarily from one protein, which we name chimallin (ChmA). Combining cryo-electron tomography of nuclear shells in bacteriophage-infected cells and cryo-electron microscopy of a minimal chimallin compartment in vitro, we show that chimallin self-assembles as a flexible sheet into closed micrometre-scale compartments. The architecture and assembly dynamics of the chimallin shell suggest mechanisms for its nucleation and growth, and its role as a scaffold for phage-encoded factors mediating macromolecular transport, cytoskeletal interactions, and viral maturation.

Breathing and tilting: mesoscale simulations illuminate influenza glycoprotein vulnerabilities.

Influenza virus has resurfaced recently from inactivity during the early stages of the COVID-19 pandemic, raising serious concerns about the nature and magnitude of future epidemics. The main antigenic targets of influenza virus are two surface glycoproteins, hemagglutinin (HA) and neuraminidase (NA). Whereas the structural and dynamical properties of both glycoproteins have been studied previously, the understanding of their plasticity in the whole-virion context is fragmented. Here, we investigate the dynamics of influenza glycoproteins in a crowded protein environment through mesoscale all-atom molecular dynamics simulations of two evolutionary-linked glycosylated influenza A whole-virion models. Our simulations reveal and kinetically characterize three main molecular motions of influenza glycoproteins: NA head tilting, HA ectodomain tilting, and HA head breathing. The flexibility of HA and NA highlights antigenically relevant conformational states, as well as facilitates the characterization of a novel monoclonal antibody, derived from human convalescent plasma, that binds to the underside of the NA head. Our work provides previously unappreciated views on the dynamics of HA and NA, advancing the understanding of their interplay and suggesting possible strategies for the design of future vaccines and antivirals against influenza.

Transmissible SARS-CoV-2 variants with resistance to clinical protease inhibitors.

Vaccines and drugs have helped reduce disease severity and blunt the spread of SARS-CoV-2. However, ongoing virus transmission, continuous evolution, and increasing selective pressures have the potential to yield viral variants capable of resisting these interventions. Here, we investigate the susceptibility of natural variants of the main protease (Mpro/3CLpro) of SARS-CoV-2 to protease inhibitors. Multiple single amino acid changes in Mpro confer resistance to nirmatrelvir (the active component of Paxlovid). An additional clinical-stage inhibitor, ensitrelvir (Xocova), shows a different resistance mutation profile. Importantly, phylogenetic analyses indicate that several of these resistant variants have pre-existed the introduction of these drugs into the human population and are capable of spreading. These results encourage the monitoring of resistance variants and the development of additional protease inhibitors and other antiviral drugs with different mechanisms of action and resistance profiles for combinatorial therapy.